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Inverting scattering experiments to obtain effective interparticle interactions is generally a 
poorly conditioned problem. L. Reatto (Phil. Mag. A 58, 37 (1986)) showed that for atomic 
liquids close to the triple point, inversions are hard because the structure closely resembles 
that of an equivalent hard-sphere fluid. Here I demonstrate that at low concentrations and for 
particles with short-ranged attractive potentials, S(k) also exhibits a very weak dependence 
on potential shape. Instead, different potentials all generate an S(k) that closely resembles 
that of the Baxter model with a similar second-virial coefficient. By contrast, in this energetic 
fluid regime, the inversion of an attractive interaction from real-space correlations such as the 
radial distribution function g(r) is well conditioned. Nevertheless, one may extract further in- 
formation from S(k) by measuring isosbestic points, values of k where the scattering intensity 
I(k) or the structure factor S(k) is invariant to changes in interaction-potential well-depth. 
These points suggest a new extended corresponding states principle for particles in solution 
based on the packing fraction, the second osmotic virial coefficient, and a new measure of 
effective potential range. 
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1. Introduction 

Extracting effective interparticle interactions from experiments is a central objec- 
tive in chemical physics because these interactions govern physical behaviour [1]. 
But like many such inverse problems, this task is complicated. Experimental data 
is not perfect: statistical fluctuations occur, and results may not be obtainable 
over a complete parameter range. Moreover, the inversion procedures are often ill 
conditioned: small differences in the initial input can result in large changes in the 
final output. 

This article concentrates on the determination of (spherically symmetric) effec- 
tive interparticle pair potentials v e f f (r) from experimental scattering data for par- 
ticles in solution. This problem is related to better studied question of how to obtain 
interatomic potentials from liquid state structure factors S(k) for simple liquids, 
where the need for very accurate experimental data, and sophisticated predictor- 
corrector inversion schemes [2] is firmly established. For simple liquids, inversions 
are subtle in part because the structure is dominated by their hard-cores [3]. For 
example, the S(k) for monotonic atomic liquids are almost quantitatively approx- 
imated by those of hard spheres (HS) [3, 4] so that the attractive interactions 
can only be determined by measurements accurate enough to distinguish the small 
deviations from HS like behaviour. 

For particles in solution, the range of relevant effective volume fractions rj is larger 
than in the better studied case of simple liquids near the triple point. In particular, 
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volume fractions are often much lower than those that lead to an entropy domi- 
nated freezing transition. At these lower volume fractions, and for short-range at- 
tractive interactions, the real-space structure can be dominated by the attractions, 
so that one can speak of energetic fluids [5]. Extraction of effective interactions from 
real-space measurements of the pair-correlation function g(r) is a well-conditioned 
problem for energetic fluids [6]. However, real-space measurements are not always 
possible or available. In this article I will show that inversions from S(k) are as 
problematic for energetic fluids as they are for entropic fluids, albeit for different 
reasons. These problems can be partially circumvented by studying isosbestic points 
- values of the wave- vector k for which S(k) is invariant under changes of the attrac- 
tive potential well depth. Isosbestic points can be used to determine the effective 
range of the v e ^(r), and also suggest an extended corresponding states principle 
for particles in solution, described by three experimentally accessible variables: the 
range defined by these points, the particle density p and the reduced second os- 
motic virial coefficient B% = B2jB^ s where B^ s is the second virial coefficient 
of a hard-sphere (HS) system. 

2. Inverting structure factors at low packing fractions 

For particles in solution, the scattering intensity I(k), measured by light, X-rays, or 
neutrons, is usually interpreted by dividing I(k) by the single particle form factor 
P(k), to obtain the structure factor S(k) = I(k)/(pP(k)) [7]; here p = N/V is the 
density. From S(k), inferences can then be made about the form of the inter-particle 
interactions v e ^{r). 

Fig. 1 helps set the stage and introduce the problem to be addressed here. Experi- 
mental data [8], taken from small-angle neutron-scattering (SANS) of a microemul- 
sion composed of small water droplets, coated with a layer of surfactant (AOT), is 
compared to some theoretical structure factors S(k). In the original letter [8] a best 
fit, using an approximate integral equation technique, yielded a hard-core diameter 
of a = QOA, a packing fraction rj = na 3 /Q = 0.075, and an attractive square well 
(SW) potential with a range A$w — 0.02cr and a well depth of /3v(a) = 3.85. The 
quality of the integral equation fit was then confirmed by independent computer 
simulations. However, as strikingly demonstrated by Fig. 1, a wide range of other 
potentials, depicted in Fig. 1(b), also lead to fits of similar accuracy. These include 
HSs with attractive SWs of two different ranges A$w, an Asakura-Oosawa (AO) 
depletion potential [9] of range Aao, an alternate simplified depletion potential of 
range Ajj showing the effects of solvation layers [10], as well as a Lennard-Jones-n 
(LJ-n) potential, defined as: 

*>-*((fr-®"). (i) 

with n = 12. The S(k) are calculated within the Percus Yevick (PY) integral 
equation closure [3, 11], and are close to the analytical PY solution of Baxter's 
model [12]. 

Improving significantly on these pioneering experiments is hard. S(k) is obtained 
by dividing I(k) through by P(k), which is only approximately determined, and 
which rapidly decays to zero for increasing ka, so that achieving better accuracies 
for a larger ka range is very difficult. As is clear from Fig. 1 (c), real-space measure- 
ments of g(r) would be much more sensitive to the effective potentials. Although 
increasingly accurate techniques being developed for the real-space measurement of 
the v e *' (r) of particles in solution [6, 13], these are often limited to certain particle 
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Figure 1. (colour online) A single set of experimental data can be interpreted by many different effective 
potentials, (a) The experimental SANS data (circles) on a microemulsion at packing fraction 77 = 0.075 is 
from [8]. The open triangles denote the S(k) from the Baxter model [12] for iJJ = —1.60, while the other 
S(k) are calculated for the potentials shown in (b). These include SWs with Agvv = 0.02 (solid lines) 
and Aj[y = 0.2 (dashed), an AO form with A^q = 0.4 (long-dashed), a generalised depletion form with 
Arj = 0.2(dottcd) and a LJ-12 (dot-dashed) potential, (c) The radial distribution function g(r) is much 
more sensitive to the potentials in this energetic fluid regime than the S(k) is. The line-styles in the three 
graphs correspond with each other. 



types or sizes. The microemulsion measured in [8], for example, would be hard to 
characterise by current real-space techniques 

The upshot of Fig. 1 is that deducing an effective potential from experimental 
scattering data at one state point is difficult; the inversion is ill-conditioned. Clearly 
more information is needed to interpret the data [14]. 

But first, what can be inferred from a single S(k)? One hint comes the Baxter 
model, which is completely determined by the packing fraction r] and the reduced 
osmotic virial coefficient I?2 or equivalently, the Baxter parameter, defined as r = 
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|/(1 — i?|). The S(k) in Fig. 1, with effective r parameters ranging from r = 0.093 
to r = 0.10 are all close to the Baxter model result with r = 0.096. To first 
order therefore, measuring a single S(k) in this regime of low packing fraction, 
commonly encountered for particles in solution, does not result in much more 
information about the effective potential than the reduced second virial coefficient 
B\. Just as found for simple liquids, where inversions are hard because the S(k) 
resemble those of HS fluids [3], here difficulties arise because they are close to the 
Baxter S(k). Distinguishing between different potentials will only be possible by 
measuring deviations from the Baxter model, and these may be very small at any 
one statepoint. 



3. Isosbestic points in structure factors 

To make further progress, measurements at other state points are necessary. This 
is done in Figs. 2(a)-(d), which depict the S(k) calculated for four LJ-n potentials 
at different temperatures. A best fit to the Baxter model S(k) at each temperature 
would result in S| as a function of temperature. But there is clearly more informa- 
tion in these curves: Each of the four v e ^ (r) results in a different set of isosbestic 
points, where the S(k) is invariant for different temperatures. Fig. 2(e) shows that 
the first isosbestic point k\o increases with increasing n, reflecting the decrease of 
the range of the LJ-n potential (1) depicted in Fig. 2(f). (Error bars in Fig 2(e) 
reflect the approximate nature of the isosbestic points.) 

These observations can be quite easily rationalised by the following simple theory: 
If the total correlation function h(r) = g(r) — 1 is split into two parts, with ho(r) = 
— 1 for r < a, and h\{r) = g{r) — 1 for r > a, then the structure factor S(k) = 
1 + ph(k) simplifies to 

S(k) = l + ^-j 1 (k) + ph 1 (k), (2) 

where ji(k) is the first spherical Bessel function, and h\(k) is the Fourier Trans- 
form (FT) of h\(r). For potentials with a hard-core, such as most of those depicted 
in Fig. 1(b), this is exact, but even for the LJ-n potentials this is a good approxi- 
mation [15]. Since ho(k) = -r^'i(fc) is independent of temperature (barring small 
effective a effects), Eq. 2 suggests that isosbestic points occur whenever h\{k) = 0. 
We note that a similar derivation was made within the MSA approximation by 
Tuinier and Vliegenthart [16], who showed that isosbestic points may also occur 
for longer ranged potentials. 

It has already been pointed out that in many experimentally relevant cases 
with low T], the g(r) are surprisingly well approximated by a simple form g(r) = 
exp(— (3v(r)) [5], as explicitly demonstrated in Fig. 3. This simple approximation, 
equivalent to taking a Mayer function [3] for h(r), works best for particles with 
short range attractive potentials. These can become quite deep, leading to large 
values of g(r) near contact, well before the system crosses a liquid- liquid or liquid- 
solid phase-line. Within this approximation, the dominant effect of varying the 
temperature is to change the amplitude of h±(r) (as demonstrated in Fig. 3(b)). To 
first order, the period of h\(k) is not affected and each ka where h\{k) = leads to 
an isosbestic point. For an infinitely narrow potential, the isosbestic points would 
= nir, but for a finite potential range there is a phase-shift, which ex- 
plains why the k n a move progressively further from mr with increasing range. The 
excellent accuracy of this simple Mayer function theory for the first isosbestic point 
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Figure 2. (colour online) (a)-(d) S(k) were calculated at r\ = 0.0576 and at different temperatures for 
each of the four LJ-n potentials v e 'f(r) given by Eq. (1). Different v e ff(r) lead to different isosbestic 
points k n cr. (e) The first isosbestic point k\a approaches tt with increasing n. A simple theory (solid line), 
described in the text, fits the data well, (f) The range of the LJ-n potentials in (a)-(d) decreases with 
increasing n. 



k\o is demonstrated by the solid line in Fig. 2(e). 

The examples depicted in Fig. 2 are at a relatively low packing fraction, but the 
isosbestic points are robust up to packing fractions of at least r\ = 0.2, as shown in 
Fig. 4(a) [17]. These points are accurately described by the simple theory, described 
above. 

Effective potentials in complex fluids may have many origins [1], so that one 
cannot simply vary the temperature to change the well depth. Nevertheless, the 
analysis above suggests that if varying a parameter leads to isosbestic points, then 
it is possible that that the well-depth of v e ^(r) is changing while the range is 
not. When combined with the variation of B2, this can help fix the form of v e $* (r). 
However, this information may still not be sufficient, since the well-depth of v e ^(r) 
may depend in an (unknown) non-linear fashion on the parameters like the temper- 
ature T, the pH, and salt or other additive concentration. Furthermore, as shown 
in Fig. 4(b), different potential shapes, picked to have the same isosbestic points, 
can generate similar S(k) for each value of B|. 
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Figure 3. (colour online) (a) The PY approximation (solid line) very accurately reproduces these molecular 
dynamics simulations [11] of a LJ-12 potential at fie = 1.6, and r) = 0.1. The simpler g(r) = exp(—f3v(r)) 
form (dashed lines) is also accurate, (b) When the temperature is changed, the dominant effect on g(r), 
calculated here with PY, is to change its amplitude. 




2 4 6 

ko 



Figure 4. (colour online) Fig (a) The position of isosbestic points barely varies with density for LJ-12, 
shown here for three temperatures f3e = 1.5, 1.33, 1.09 at each density r\ = 0.05, 0.1, 0.15, 0.2. The line with 
the arrow denotes k\a for increasing r). (b) Three potentials leading to the same isosbestic points: SW 
(A sw = 0.26) (solid lines), AO (A AO = 0.6) (dashed lines) and LJ-22 (dot-dashed lines) for r) = 0.05 and 
f?| = 0, —0.719, —1 (plots in descending order of B2 at k = 0). 



4. An extended corresponding states principle 



This similarity of the S(k) ties in with the work of Noro and Frenkel (NF), who 
proposed that many properties of particles in solution could be understood from 
an extended corresponding state principle based on the variables rj, B^-, and an 
effective well depth j3e [18]. However, the three potentials shown in Fig. 4(b) do 
not have the same well-depth, while still showing similar S(k). This suggests an 
alternative corresponding states principle based on the variables rj, B%, and an 
effective range A e jj related to the isosbestic points. In fact, as NF point out, the 
potential range is not always uniquely defined when comparing different v e ^ (r)'s. 
They suggested a non-linear mapping based on B^ to derive an effective SW range 
for each potential they consider. 

Here I define the effective range A e ff of a given potential as being that of an 
equivalent SW fluid with the same same k±a. A simple and accurate approximation 
of k\o for the SW fluid at low concentrations is given by k\o ~ 7r/(l + A e jf/2 + 
Agj^/12) which follows from an expansion of the exact solution to the simple theory 
of isosbestic points. Thus the value of the effective range is then well approximated 
by: 

A ef f = 24 ( - — V - 12. (3) 



By comparing the k\a values for different potentials, some simple empirical rules 
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Figure 5. Defining effective range parameters from isosbestic points. For each potential type a simple 
mapping described in the text was used to derive an effective square well range A e ff- For each potential, 
a PY calculation was used to find the isosbestic points, and these are shown in the plot for four different 
kinds of potentials. The long-dashed line denotes an accurate analytical approximation for the isosbestic 
points of square- well fluids. If one does not know what the potential is, then ^ e ff can De directly estimated 
from the k\a measured in experiments. 

can be derived: For the LJ-n potentials this criterion reduces to the distance r — a 
where (3v{r) reaches 1/4 of its maximum depth, for the AO potential the effective 
SW range is ^ e ff — 0.42A J 4o, and for the hard-core Yukawa potential, with an at- 
tractive tail of the form v(r) = eexp(— r/A) for r > a, the mapping is A e ff — 1-15A. 
The results are depicted in Fig. 5, and show that these simple linear mappings suc- 
cessfully define an effective range for each potential type. 

Another argument for including the A e jj (instead of the well-depth) as one of 
the relevant 3 variables is that the range determines the topology of the phase- 
diagram [18]. The critical range below which the fluid-fluid transition becomes 
metastable to the fluid-solid line is at A e jy « 0.15 for all four potentials (SW, AO, 
Yukawa and LJ-n), something consistent with what NF found. In addition, A e jj 
is directly experimentally accessible through k\a. For example, taking advantage 
of the proximity of a metastable fluid-fluid critical point to enhance nucleation, a 
proposal made for protein solutions [19], would involve tailoring solution conditions 
such that the measured isosbestic point is just above k\o = 2.92, so that A e j^ is 
just below 0.15. 



5. Discussion 

In summary then, extracting effective potentials v e ^(r) from S(k) at the lower 
packing fractions typically encountered for particles in solution is a poorly condi- 
tioned problem. Whereas in simple liquids the HS model describes the dominant 
features of the S(k), here the Baxter model appears to be the fundamental un- 
derlying model for S(k) around which different attractive potentials only induce 
a mild perturbation. More information can be obtained by studying the isosbestic 
points, which help determine the range A e ^ of the potentials. Together with rj and 
I?!, A e jj can be used to define an extended corresponding states principle. Parti- 
cles in solution with similar values of these three parameters should have similar 
properties, such as the relative stability of the fluid-fluid and fluid-solid binodals. In 
other words, many solution properties could be deduced without having to actually 
invert to an explicit form of v e ^ (r). 
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Nevertheless, It should always be kept in mind that the effective pair potential 
that reproduces the structure, even though it is unique [20], may not be the same 
as the effective potential that reproduces certain thermodynamic properties. Such 
representability problems [21] are typically important when there are large three- 
body interactions, or when the underlying system has anisotropic interactions [22]. 
For systems with truly short-ranged potentials, one doesn't expect strong many- 
body effects. Nevertheless, anisotropic interactions may be common for colloidal 
suspensions. They lead to important deviations from spherically symmetric poten- 
tials [23-26], could lead to representability problems if a potential is inferred from 
inverting structure factors. 

In addition, particles in solution are often polydisperse, an effect that has been 
ignored in the current analysis. Polydispersity can manifest in multiple properties, 
i.e. in the effective diameters as well as in the interaction strengths. Clearly this 
adds a different dimension of complexity, and it would be interesting to see how 
isosbestic points would be affected by polydispersity. Encouragingly, structure fac- 
tors have been solved for polydispersity in systems like the Baxter model [27], so 
that progress may be expected on this front in the future. 

Finally, numerous examples of isosbestic points in S{k) can be found in the liter- 
ature, ranging from colloids with short-range sticky coats [28], to colloid-polymer 
mixtures [29], to globular proteins [30], to magnetic colloids [31] and to colloid- 
micelle mixtures [32]. Interestingly, these points should also appear in I(k). 
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